home *** CD-ROM | disk | FTP | other *** search
/ IRIX Base Documentation 2001 May / SGI IRIX Base Documentation 2001 May.iso / usr / share / catman / p_man / cat3 / complib / dsbev.z / dsbev
Encoding:
Text File  |  1998-10-30  |  3.7 KB  |  133 lines
  1.  
  2.  
  3.  
  4. DDDDSSSSBBBBEEEEVVVV((((3333FFFF))))                                                            DDDDSSSSBBBBEEEEVVVV((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      DSBEV - compute all the eigenvalues and, optionally, eigenvectors of a
  10.      real symmetric band matrix A
  11.  
  12. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  13.      SUBROUTINE DSBEV( JOBZ, UPLO, N, KD, AB, LDAB, W, Z, LDZ, WORK, INFO )
  14.  
  15.          CHARACTER     JOBZ, UPLO
  16.  
  17.          INTEGER       INFO, KD, LDAB, LDZ, N
  18.  
  19.          DOUBLE        PRECISION AB( LDAB, * ), W( * ), WORK( * ), Z( LDZ, * )
  20.  
  21. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  22.      DSBEV computes all the eigenvalues and, optionally, eigenvectors of a
  23.      real symmetric band matrix A.
  24.  
  25.  
  26. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  27.      JOBZ    (input) CHARACTER*1
  28.              = 'N':  Compute eigenvalues only;
  29.              = 'V':  Compute eigenvalues and eigenvectors.
  30.  
  31.      UPLO    (input) CHARACTER*1
  32.              = 'U':  Upper triangle of A is stored;
  33.              = 'L':  Lower triangle of A is stored.
  34.  
  35.      N       (input) INTEGER
  36.              The order of the matrix A.  N >= 0.
  37.  
  38.      KD      (input) INTEGER
  39.              The number of superdiagonals of the matrix A if UPLO = 'U', or
  40.              the number of subdiagonals if UPLO = 'L'.  KD >= 0.
  41.  
  42.      AB      (input/output) DOUBLE PRECISION array, dimension (LDAB, N)
  43.              On entry, the upper or lower triangle of the symmetric band
  44.              matrix A, stored in the first KD+1 rows of the array.  The j-th
  45.              column of A is stored in the j-th column of the array AB as
  46.              follows:  if UPLO = 'U', AB(kd+1+i-j,j) = A(i,j) for max(1,j-
  47.              kd)<=i<=j; if UPLO = 'L', AB(1+i-j,j)    = A(i,j) for
  48.              j<=i<=min(n,j+kd).
  49.  
  50.              On exit, AB is overwritten by values generated during the
  51.              reduction to tridiagonal form.  If UPLO = 'U', the first
  52.              superdiagonal and the diagonal of the tridiagonal matrix T are
  53.              returned in rows KD and KD+1 of AB, and if UPLO = 'L', the
  54.              diagonal and first subdiagonal of T are returned in the first two
  55.              rows of AB.
  56.  
  57.      LDAB    (input) INTEGER
  58.              The leading dimension of the array AB.  LDAB >= KD + 1.
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. DDDDSSSSBBBBEEEEVVVV((((3333FFFF))))                                                            DDDDSSSSBBBBEEEEVVVV((((3333FFFF))))
  71.  
  72.  
  73.  
  74.      W       (output) DOUBLE PRECISION array, dimension (N)
  75.              If INFO = 0, the eigenvalues in ascending order.
  76.  
  77.      Z       (output) DOUBLE PRECISION array, dimension (LDZ, N)
  78.              If JOBZ = 'V', then if INFO = 0, Z contains the orthonormal
  79.              eigenvectors of the matrix A, with the i-th column of Z holding
  80.              the eigenvector associated with W(i).  If JOBZ = 'N', then Z is
  81.              not referenced.
  82.  
  83.      LDZ     (input) INTEGER
  84.              The leading dimension of the array Z.  LDZ >= 1, and if JOBZ =
  85.              'V', LDZ >= max(1,N).
  86.  
  87.      WORK    (workspace) DOUBLE PRECISION array, dimension (max(1,3*N-2))
  88.  
  89.      INFO    (output) INTEGER
  90.              = 0:  successful exit
  91.              < 0:  if INFO = -i, the i-th argument had an illegal value
  92.              > 0:  if INFO = i, the algorithm failed to converge; i off-
  93.              diagonal elements of an intermediate tridiagonal form did not
  94.              converge to zero.
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.